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Abstract 

We give the results for all the one-loop propagators, including finite parts, in the Coulomb 
gauge. In finite parts we find new non-rational functions in addition to the single logarithms 
of the Feynman gauge. Of course, the two gauges must agree for any gauge invariant 
function. 
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1. Introduction 

The non-covariant axial and Coulomb gauges have more direct physical interpretation 
than the covariant gauges because their propagators are closely related to the polarization 
states of real spin-1 particles. The relevant diagrams in the Coulomb gauge are not plagued 
by ghosts. Also the time-time component of the gluon propagator provides a long-range 
confining force [1], [2]. The Hamiltonian for non-Abelian gauge theory in the Coulomb 
gauge has been known for some time in its continuum version [3]. The Coulomb gauge 
in the Hamiltonian formalism is manifestly unitary. The main point in its favour is that 
problems concerned with the definition of the axial gauge integrals like 

do not appear in the definition of integrals like 



/ 



d 4 k^-.. (2) 



in the Coulomb gauge. However, there are disadvantages. The naive Coulomb gauge 
Feynman rules in non- Abelian gauge theory give rise to ambiguous integrals, in addition 
to the usual ultra-violet divergences [4]. At one loop order and above there are integrals 
like 

r d 3 P r d Po po i 

J (2-KfJ (2tt) pi - P 2 + i V (P-K)*' [ ] 

There is no regularization procedure for the energy divergence in po within the standard 
dimensional regularization scheme. This integral and similar more complicated divergences 
in higher order diagrams have been the subject of the study [5] , [6] , where systematic can- 
cellations have been found. However, no general proof exists that controls all divergences 
[7]. Formally such integrals are assigned value zero. The Coulomb gauge has been exten- 
sively studied in the phase space formalism by D. Zwanziger [8] in the Euclidean space. 
The ultra-violet divergent parts of the proper two-point functions have been calculated and 
found to observe the Ward identities. In addition, a more powerful Ward identity holds in 
the Coulomb gauge than is available in covariant gauges. In this paper we give the results 
for the complete propagator to order g 2 including finite parts in Minkowski space. 
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2. The Coulomb gauge in the phase-space formalism 

We use the phase-space formalism in order to avoid the ambiguous integrals like (3). Let 
the generating functional of the Green's functions be 

Z(j, J) = J d[f]d[A][J^ + j^f^} expH J d 4 xL], (4) 

where J, j are sources, L the Lagrangian density and / and A are the fields [9]. 

L = -\f;j a " v + \n v F^ (5) 

and 

F; u = d^Al - d v Al - gf ahc AlA c u . (6) 

The Greek indices run from to 3 and Latin indices denote spatial dimensions (i = 1, 2, 3), 
a, 6, c, are colour indices. Instead of setting the source j to zero as in Lagrangian formalism, 
we keep some components of j. We write out L as 

L = -\{ft 3 ? + \{f^f + \ft 3 F* - (7) 

and set jij = in order to perform the Feynman integral over We denote /oi = Ei. 
The Lagrangian becomes 

L = -\(F? j ) 2 + \{E1) 2 -E?F$ l (8) 

where 

W = E t[doA1 - d t A a - gt hc A\A^. (9) 

The field Ei is the momentum conjugate to Ai. We have L expressed in terms of momenta 
and linear terms in derivatives. It is first order in time derivatives. Adding the gauge 
fixing term 

^(^) 2 (10) 

we can deduce the propagators from the quadratic part of the Lagrangian. The propagator 
matrix is the inverse of the matrix of these quadratic parts, and is given by the following 
7x7 matrix, whose rows and columns are labeled by A\, A 2 , A3 ; A ; E±, E 2 , E 3 : 
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A 



Ai -T l3 /k 2 + aL lJ /K 2 


ak Kj(K 2 ) 2 


-ik T in /k 2 


A ak Kj/(K 2 y 


1/K 2 + ak£/(K 2 ) 2 


iK n /K 2 


E m —ikoTmj/k 2 


iK m /K 2 


-Tmr,K 2 lk 2 


where 






Tij = Sij Lij , 


= KiKj/K 2 , 




k 2 


= k 2 - K 2 . 


(11) 


The Coulomb gauge propagators are obtained by setting a = 0. 




3. The proper two-point functions 






The method of evaluation of Coulomb g 


;auge integrals is explained 


in Appendix A. Here 


we list the results. The constants used throughout the paper are 





e = 4-d (12) 
where d is the dimension of space-time and the coupling parameter is 

• 2 

cSab = -r^CcSab (13) 



The transverse gluon two-point function 

There are two non- vanishing graphs contributing to the transverse gluon propagator. The 
graph shown in Fig.l gives 

^=cT(l)(5r f x2-.(l + ^) 

xf^A' 2 ^ - ^K,Kj + j^cKiKj] (14) 
The graph in Fig. 2 contributes 

Tz iAj = c(MKiKj + NK 2 5 lJ ) (15) 

where 
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31 4, ,-k 2 -ir>^ 3 , if 2 

1 K A h 2 h 2 h 2 h A 

+ 4 [ ~ fcg - + 18 ^T + ^ + K ] ln fc - if + Tr, X ln (-fc 2 - in) 
L -[3^-5^ + 9^-^] In- 



2i^ 2L if 2 fc 2J {-k 2 -ir,) 

In 9 t 4 1 t 4 ?T 4 



K 2l "K 2 ' 5" ' "if 2 fc 2 

+ 22^|-16- ---+4x — (16) 
if 2 9 15 225 V ; 



4 L fc 3 1 if 3 /c 2 if 2 ^ J k -K + iri (-k 2 -i7]) 

I ... . A- 1 A ' if 2 



ln2,_. .o /c 4 if 4 11 



fr 2 



3^- 


K 4 


if 2 


k 2 


A; 4 


K 4 


3 if 2 " 




+ 10 + 


2 
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■jr- ' ™ ' 9 (17) 

The non-rational structure D which appears in the results for the proper two-point 
functions is in the integral form 



D 



= l!< h kZ-x(K>-ir l ) H1 - X) - (18) 
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In the region k > K 



1 ko-K + ir] ko + K-ir] 

k K X ^ko + K-irj' ^ko-K + ir]' 

ko + K-ir] k 2 + ir] k + K - ir] 

+ ln x In — - — Z7rln -\. (19) 

ko-K + ir] K 2 ko-K + ir]' 



In the region K > ko 



1 K -kp-ir] K + kg — ir] 

d = Tn7\ Ll ^-[m. — ~ ~ Ll2 ^i? — i — ~ 

koK K + ko — ir] K — ko — ir] 



K + ko - ir] ^ -k 2 - ir] ^ ^ K + k - ir] - 
K - k - it] ko K - k -ir]~ 



* * (20) 



kpK K — ir] K — ir] kpK (—k 2 — irj) 



where 



Li 2 (x) = - r in(1 Z) dz (21) 
Jo z 

is the Spence function and ko and K are the lengths of the respective vectors. The two 
expressions for D in (19) and (20) are connected as analytic continuations of each other 
with the relation (Bl). 

The AiAo transition 

The whole contribution to the AiAo transition to order g 2 comes from the graph in Fig. 3a. 



T AiA ° = ck Ki x Z (22) 



+^r(2^ + fc 2 ) x D 



k 2 /^,o ,n Sl k + K-iri , K 2 
■(2k 2 + k 2 ) In -f ^ x In 



2/c if 3 ko-K + ir] (-k 2 -iri) 

„ k _ t£ k _ _ kr\ 53 , x 

-V n R^- 6 ^ ln2+6 7^-y < 23 > 

The graph shown in Fig. 3b contains integrals like the one in (3). Formally such integrals 
are assigned value zero. 
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The time-time component of the two-point function 

Two graphs contribute to the A A function. The graph shown in Fig. 4a gives 

= c { r ( i)(^)-» x (\ k l + \# + JL*» + l*g + ^) 



2 



2 



+ lk 4 x D 



fc 4 , ko + K-ir] , if 2 
In x In ■ 



2/c if k - K + ir] (-k 2 - irf) 

~ fe ° 2in (-X^ ~ 2(in2 ~ i)/g ° 2} (24) 



The graph in Fig. 4b contributes 



+/c 2 /c 2 x L> 

fco/c 2 , k + K-ir] , 
In x In 



K ko — K + ir] (—k 2 — irj) 
K 2 

-(2k 2 + K 2 ) In — - 2{2k 2 + K 2 )(ln2 - 1)} (25) 

(— k A — ir\) 

We can verify that the complete proper two-point functions satisfy the 't Hooft identity 
[10] 

klT AoAo - 2k KiT AiAo + KiKjT AiAj = (26) 
and the stronger Zwanziger identity [8] 

k T AoAo = KiT AiA °. (27) 

The remaining graphs contain the conjugate field Ei as the external leg. 

EiAj graph 

The graph in Fig. 5 vanishes as the energy divergence. 
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EiA graph 

The graph in Fig. 6 gives 



F ^A = c( _)-5{_ _ 2"T(-)(- + -p)} x (2z*Q) (28) 



L 3 x 2 /x 3 9 



EiEj graph 

The graph in Fig. 7 contributes 



r E iEj = _ 2c (^)-§ x {2" e (- + — )r(-)^- - i^El} (29) 
V 2 3 9 2 3 K 2 y J 



4. The gluon propagator to order g 2 

We form the 7x7 matrix of free and order g 2 proper two-point functions. The inverse of 
this matrix gives the gluon propagator to order g 2 . 

The AqAo propagator 

The time-time component of the gluon propagator to order g 2 is [11] 

K K 



or explicitly 



D A A = , K 2y2 { ^L V (-)K 2 - -K 2 In ( ^ 9 ^ - 2K 2 In % 

1 3 V 3 ^ 2 ^ 2 



+\k 2 {k 2 + 2k%) x £) 
:(A; 2 + 2/c 2 ) In — — i x In ■ 



2fc tf u ' ko-K + i V (-k 2 -i V ) 

-(3k 2 - K 2 ) In - ,f 2 , - - (6A; 2 + 2K 2 ) In 2 + 6k 2 + ^K 2 } (31) 
{—k A —ir]) 9 

The ultraviolet divergent part of (31) gives the gauge invariant Coulomb field renormal- 
ization factor [12]. 

The AiAj propagator 

The transverse gluon propagator to order g 2 is 

1 k 2 

TlAiAj rp -p Am A n rp 0_rp -pE m E n rp f'QO'i 

J - / — ^- L am 1 - - L nj ^- L am 1 - - L nj V 3 ") 
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or explicitly 



c KiKj. r _,e, 4 K 2 1 ,-k 2 -ir], 

k A +17] K z 2 6 fX 2 6 (I z 



^r ^ 2 + fe o + ^(14 + 3^)1 In fco + K ~^ x In ^ 



1 A; 2 if 2 if 2 



/c 2 if 2 37 N1 92, , , 

! -(3^ + — + y) ln2 - y> ( 33 ) 



5. The Slavnov- Taylor identity 

Although ghosts are absent from the S-matrix elements they are necessary to formulate 
the Slavnov- Taylor identities [13], [14]. Diagramatically they are shown for the self-energy 
in Fig. 8. Algebraically they are 

k T AoA i - KiT AiA J = (K 2 5ij - KiKj)T CAi (34) 

The diagrams involving ghost-source vertices on the right-hand side are shown in Fig. 9a 
and Fig. 9b. The diagram in Fig. 9a vanishes as the energy divergence in p . The diagram 
in Fig. 9b contributes 

r CA > = -2c(^)-iK l{ - 4 - + 2"T(i)(| + ^)} (35) 

so the identity is satisfied trivially as implied by (26) and (27). 
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6. Discussion 

We have checked the consistency of the Coulomb gauge to order g 2 including finite parts. 
The time-time component of the gluon propagator in the Coulomb gauge is believed to 
provide a long-range confining force. There are two interesting limits of eq.(31). In the 
Zwanziger picture [8] g 2 D 00 gives the instantaneous part Vz(R), which is called the color- 
Coulomb potential. (Here D 00 is the time-time component of the gluon propagator). The 
instantaneous color-Coulomb potential Vz(R) at large R may serve as an order parameter. 

Kcoui = lim (36) 

A non-zero value of Kcoui would be the signal for color confinement. The potential is 
separated out in momentum space by 

V Z {K) = lim g 2 D A ° A °(k ,K) (37) 

fco— >oo 

where we have written Vz(K) for the Fourier transform of Vz(R). The limit ko — > oo of 
eq.(31) is 

lim ^o (fc0)K) = ^ r( | ) _^ ln ^!_, 7r _| ln2+ l^_ 21n ^ } (38) 
fc ->-oo K z 3 2 3 \i A 3 9 Kq 

and it is not independent of ko. There appears to be a difference in the dominant term 
in (38) and equation (37) of A. Cucchieri and D. Zwanziger [15]. This difference arises 
because of the statement near the end of Appendix B in [15] that h is finite, and "as 
a result" I2 vanishes in the limit ko — > 00. However, the finitness of I2 does not imply 
anything about the behaviour as ko — > 00. In fact, on calculating I2, we find that the 
dominant term as ko — > 00 is -4/3ln(k 2 /K 2 ). With this value, there is no contradiction 
with (38) in this paper. 

Although the limit as k — > 00 is not finite, A. Cucchieri and D. Zwanziger [15] have 
argued that an unambiguous instantaneous part may be defined by using renormalization 
group arguments. 

The limit ko — > is naturally related to the definition of the quark-antiquark potential. 
It follows from considering a rectangular Wilson loop with sides of length T in the time 
direction (where T — > 00) and L in the space direction. In the Coulomb gauge the main 
contribution comes from the D o component of the propagator (where ko — > 0) attached 
to the two time-like sides. The k — > limit of eq.(31) is 

&^<^ = jMr r <5>-T to ? + T> (39) 

leading to the quark-antiquark potential 

V(R) = -2n 2 g 2 M±{l + + y 7 + y ln(^) 2 ]} (40) 

where 7 is the Euler's constant, g r (lj) is the running coupling constant. If we assume the 
relation 

R x n = 1 (41) 
<7 r (//) becomes R dependent. We suppose that the exact <7r(j=>) tends to zero as R — > 
and g r {ji) ~^ 00 f° r R °°- 
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Appendix A 

We use two basic integrals for evaluation of the Coulomb gauge integrals. 



p 2 + if] (k — p) 2 + if] 
Un J dy J d 3 ~ e P{P 2 -2P-Ky- yk 2 + y 2 kl - it/} - ' (Al) 



B = I d 4 -p- P0 



(p 2 + irj)[(k - p) 2 + irj\ 
= ^ink J ydy J d 3 ~ e P{P 2 - 2P ■ Ky - yk 2 + y 2 k 2 - it/}"* (A2) 
As an example we evaluate the integral 

Xa = [ d*-« p -25_ - L ^ (^3) 

3 J F p 2 + ir] (k-p) 2 + ir] P 2 v ; 

Applying (A2) 

1 f 1 f P P 1 

Xa = -ink Q / ydy / d 3 ~ e P^- T (AA) 

3 2 U 7o 7 P 2 [P 2 -2P-Ky-yk 2 + y 2 k 2 -n 1 ] 3 2 V ' 

Combining the denominators with the Feynman parameter x 
X l3 = iko^T(-) J d.r.r- J ydy J d 3 ~ e P 



[P 2 — 2P ■ Kxy — xyk 2 + y 2 xkg — irjx] 2 

(A5) 

Now it is easy to perform the d 3 e P and integration over the parameter y giving Xij in 
the integral form. 

C- X X l3 = \V{\){K 2 )-H tj - + \6iA™ ~ In 2) 

fix ^± \h 2 S f\ x ~ 2 2 f 1 a* K 2 {l-x) 

K 2 > k iJ Q dx k 2 -x(K 2 -i V ) + J Q dX {k 2 -x(K 2 -i V )} 2 (-k 2 -i V ) } 

KiKj r0r l f 1 , x^ , o f 1 , x^ 



- kl{ 2] dX k 2 -x(K 2 - l r 1 ) +k2 j dX \ 



K 2 01 2J k 2 -x(K 2 -ir]) J [k 2 - x{K 2 - irj)]- 
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, ,4 f\ ^ (l-x)K\ 

+fc Jo [** - X {k* - i V )]* m - ir,y 

where 

C = ik Q TT^. (A6) 

The integrals in (A6) are 



1 , 2_ _fcn kp + K -irj 

k 2 - x(K 2 - i V ) K 2 K 3 ko — K + irj { ' 

1 , ^ _J _2 i k + K -ir] 

n i - x{K 2 - irj)} 2 k 2 K 2 2k K 3 U k -K + i V 

^ 1 1 1 ] k + K-i V 

Jo [k 2 - x(K 2 - irj)] 3 2K 2 k* Ak 2 K 2 k 2 8k 3 K 3 k - K + irj { } 



o 



-/ 

Jo 



x-2 , 2 f x^ K (1 — x) 

dX— o -—-^ ; — r + k I dX-r-7; = — TTTT hi 



k,Q — x(K 2 — irj) J [k 2 — x(K 2 — it])] 2 (—k 2 — irj) 



2 „ „ ,v r 1 k 2 ko + K-tri, ^ K 2 k 2 „ 

= — (ln2-l)+— 7 ln — xln- 5 r (A10) 

K 2V ; L K 2 2£; if 3 ko-K + irj 1 (-k 2 - irj) 2K 2 y 1 

/c 2 + if 2 1 1 k + K -irj 1 

~ 2£; 2 iW ~ 2iW ~ 2Xfc fc 4 A; - K + zr? ~ 8K 2 ! 2 ( } 

where D was defined in (18) and the explicit result given in (19) and (20). 
Appendix B 

The expressions for D in (19) for ko > K and in (20) for K > ko ought to be connected by 
analytic continuation. It is easy to see this happens using the following relation between 
the Spence functions. 

r . ,x — l + irj, T . ,x + 1 - irj. T . ,1 + x — ir). T . ,1 — x — ir). 

Li 2 {^— - Li 2 { 7—— + Li2{- - Li 2 {— -) 

x + 1 — irj x — 1 + irj 1 — x — irj 1 + x — irj 

+2Li 2 (—x — irj) — 2Li 2 (x + irj) + In x ln(x 2 ) 

x — 1 + irj 

-in In X + 1 - 171 + in \n(x 2 ) + tt 2 = (51) 
x-l + irj v ' v ' 



where 



Lr 2 (x) = - J* ln{1 z Z) dz (B2) 
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Figure Captions 

Fig.l. The transverse gluon self-energy graph. The dashed line is the transverse gluon A i: 
the dotted line represents the instantaneous Coulomb field A and the continuous line is 
the Ei field conjugate to the transverse field A t . 

Fig. 2. The transverse gluon self-energy graph. The dashed line is the transverse gluon 
field Ai. 

Fig. 3a. The AiA two-point function. The dotted line is the instantaneous Coulomb field 
Aq, the dashed line represents the transverse field Ai and the solid line is the conjugate 
field Ei. 

Fig. 3b. The AiA two-point function. The graph is suppressed as an energy divergence. 

Fig. 4a. The time-time component of the gluon self-energy to order g 2 . The dotted lines 
represent the instantaneous Coulomb field A . The continuous line is the Ei field conjugate 
to the transverse field A t . The propagators inside the loop are the EiAj transitions specific 
to the Coulomb gauge. 

Fig. 4b. Self-energy graph to order g 2 . The dotted line is the A field, the dashed line is 
the transverse propagator and the solid line is the EiEj propagator. 

Fig. 5. The transition between the transverse gluon field and its conjugate field Ei. 
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Fig. 6. The transition between the Coulomb field A and the conjugate field Ei. 
Fig. 7. The conjugate field self-energy. 

Fig. 8. The Slavnov- Taylor identity for self-energy graphs. The wavy lines stand for Yang- 
Mills particles and double lines for ghosts. The symbol on the left wavy line stands for 
the replacement of a polarization vector e /x (fc) by fc M and k 2 need not be zero. The cross 
denotes the action of the tensor {k^k v — /c 2 5 Mi ,). The circle represents the set of all relevant 
Feynman graphs. 

Fig. 9a. Diagram with an open ghost line. The source v n of the E m field has the vertex 
gf ahc E h n C c v n . The ghost propagator is and it is represented with the double line. 

Fig. 9b. Diagram with an open ghost line. The source uf of the transverse gluon field has 
the vertex gf abc 5ij. 
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